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INTRODUCTION 


A  variety  of  calculational ,  analytical  and  literature  survey 
studies,  results  of  which  are  reported  here,  were  undertaken  by  Atmospheric 
Science  Associates  to  assist  the  Cloud  Physics  Branch,  Meteorology  Division, 
Air  Force  Geophysics  Laboratory  in  the  design  of  a  cloud  water  content  meter 
(WCM)  and  a  hydrometeor  spectrometer  (HS).  Both  of  these  instruments  are 
to  be  airborne  on  cloud  physics  research  airplanes.  The  critical  measure¬ 
ment  element  of  both  is  a  flat,  circular  annular  surface  which  is  rotating 
at  high  speed  and  is  oriented  normal  to  the  airplane  flight  path.  Thus, 
hydrometeor  impact  trajectories  are  approximately  normal  to  the  annular 
surfaces.  Both  instruments  use  momentum  effects  of  collected  hydrometeors 
on  rotation  speed  to  measure,  in  the  case  of  the  WCM,  cloud  water  con¬ 
tent,  and  in  the  case  of  the  HS,  masses  of  individual  hydrometeors.  Rotation 
of  the  annuli  is  caused  by  flow  of  air  against  vanes,  or  torque  blades  as 
they  are  called  here,  which  are  mounted  on  the  annuli  surfaces.  Figure  1 
is  an  artist's  sketch  of  the  geometry  of  the  WCM. 

Our  studies  have  been  of  three  types: 

1.  Flow  field  and  water  drop  trajectory  calculations. 

2.  Theoretical  analysis  of  torque  on  the  torque  blades. 

3.  Hydrometeor  impact  effects. 

Details  are  presented  in  the  three  chapters  which  follow.  These  studies 
are  but  a  fraction  of  the  total  design  effort,  the  bulk  of  the  work  having 
been  done  in-house  at  the  Air  Force  Geophysics  Laboratory.  The  part  of 
Atmospheric  Science  Associates  in  the  overall  effort  was  to  assist  in  the 
solution  of  various  technical  problems,  or  answer  various  technical  questions, 
as  these  arose  during  the  design  process. 

Most  of  the  effort  has  been  devoted  to  the  WCM,  and  the  design  of 
this  instrument  has  progressed  to  a  final  stage.  There  has  been  time  for 
only  preliminary  studies  on  the  HS. 
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Figure  1.  Artists  rendition  of  the  water  content 
meter.  (Not  to  proportional  scale.) 
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FLOW  FIELD  AND  WATER  DROP  TRAJECTORY  CALCULATIONS 


GENERAL 


Flow  field  and  trajectory  calculations  were  performed  by  use  of 

the  codes  developed  by  Norment^).  The  codes  calculate  potential  flow  in 

three  dimensions  about  bodies  of  arbitrary  shape  by  use  of  the  method  of 
(21 

Hess  and  Smith'  ' .  Trajectories  of  water  drops  of  any  size  are  calculated 
through  the  flow  fields  by  numerical  integration  of  the  water  drop  equaticr 
of  motion  in  which  drag  relations  derived  from  experimental  data  are  used. 
Details  are  given  in  reference  1. 

Since  potential  flow  is  calculated,  effects  of  friction  are  not 
included,  and  the  existence  of  a  boundary  layer  is  simply  ignored.  While 
details  of  the  boundary  layer  flow  would  be  necessary  for  accurate  calcula¬ 
tion  of  dynamic  pressure  forces  near  the  surface  of  the  annuli  and  blades, 
the  effort  required  would  have  been  beyond  the  scope  of  this  program.  More¬ 
over,  most  of  the  critical  information  required  for  the  instrument  design  can 
be  obtained  without  considering  the  details  of  the  boundary  layer  flow. 

To  calculate  flow  fields  about  the  various  bodies  of  interest,  it 
is  necessary  to  develop  digital  descriptions  of  the  surfaces  of  the  bodies 
in  terms  of  appropriate  configurations  of  many  small,  contiguous,  plane 
panels  or  elements.  Computer-drawn  plots  of  the  element  configurations 
are  given  below,  along  with  the  results  of  the  flow  field  and  trajectory 
calculations. 

Since  friction  is  ignored,  rotation  of  the  annuli  also  is  ignored. 
This  being  the  case,  accounting  for  flow  about  the  rotating  blades  is  some¬ 
what  of  a  problem.  In  some  of  the  calculations  discussed  below  the  blades 
are  ignored,  but  for  the  final  configuration  of  the  WCM  they  are  included 
in  approximate  form. 


1.  H.  G.  Norment,  "Calculation  of  Water  Drop  Trajectories  to  and  About 
Arbitrary  Three-Dimensional  Bodies  in  Potential  Airflow,"  Atmospheric 
Science  Associates,  NASA  CR  3291  (August  1980).  N80-28302 

2.  J.  L.  Hess  and  A.  M.  0.  Smith,  "Calculation  of  Non-Lifting  Potential 
Flow  About  Arbitrary  Three-Dimensional  Bodies,"  McDonnell  Douglas 
Report  E.  S.  40622  (15  March  1962).  AD-282  255 
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SIMPLE  DISC 


Figure  2  shows  the  digital  description  used  to  represent  a 

simple  disc  collector  of  radius  R  =  3  cm.  In  this  case,  as  in  those  to 

follow,  the  downstream  housing  is  represented  by  a  crudely  resolved 

cylinder,  while  the  collection  face  is  sufficiently  well  resolved  to  allow 

for  determination  of  the  flow  field  near  its  surface. 

Water  drop  trajectories  are  shown  in  Figures  3.  Figure  3a  can 

be  considered  to  show  flow  streamlines.  Impaction  efficiencies,  E,  cal- 

★ 

culated  from  these  trajectories  are  shown  on  the  figures  .  Figure  4 
shows  the  flow  field  near  the  surface. 

DISC  WITH  PARABOLIC  HUB 

The  simple  disc  collector  was  modified  by  addition  of  a  non¬ 
rotating,  cylindrical ly  symnetric  parabolic  hub,  which  at  its  extreme 
upwind  extension  ends  in  a  needle-like  point.  The  equation  of  the  parabola 
is 

(r/r2  -  5/6)2  =  -  l(x/r2  -  1/12)  (1) 

where  r  is  radial  distance  from  the  disc  axis,  x  is  distance  in  the  direction 
of  the  free  stream  flow,  with  origin  at  the  disc  face,  and  r2  is  the  disc 
radius  (3  cm). 

Figure  5  shows  the  crossection  plane  through  the  device  that  in¬ 
cludes  the  symmetry  axis.  The  hub  Intersects  the  face  'f  the  disc  at 
coordinates  (x/r2  =  0,  r/r2  =  2/3),  and  the  focus  of  the  parabola  is  at 
point  (x/r2  *  0,  r/r2  =  5/6).  The  annular  section  of  the  disc  surface  be¬ 
tween  r/r2  =  2/3  and  r/r2  =  1  Is  the  collection  (i.e.,  measurement)  surface, 
and  is  rotating  at  high  speed. 

^Impaction  efficiency  is  the  ratio  of  particle  flux  actually  incident  on 
the  collection  surface  to  particle  flux  in  the  free  stream. 


Figure  3b.  Trajectories  of  water  drops  of  radius  5  mn  to  the 
disc  collector.  E  Is  impaction  efficiency. 


Figure  3c.  Trajectories  of  water  drops  of  radius  T5  um 
to  the  disc  collector. 
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Figure  3d.  Trajectories  of  water  drops  of  radius  35  um 
to  the  disc  collector. 
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Figure  4.  Flow  field  near  the  surface  of  the  disc  collector 
shown  In  Fig.  2. 
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Figure  5.  Crossectlon  profile  of  the  disc  collector 
with  parabolic  hub.  The  x  axis  Is  a  rotation  symmetry 
axis. 


Figure  6  shows  the  digital  description  used  to  represent  this 
version  of  the  WCM.  Figures  7  show  water  drop  trajectories,  with  those 

in  Figure  7a  being  essentially  coincident  with  flow  stream  lines.  Inpaction 
efficiencies,  E,  calculated  from  the  trajectories  are  shown  on  the  figures. 
Figure  8  shows  the  flow  field  close  to  the  surface  of  the  rotating  annulus, 
and  Figure  9  shows  contours  of  constant  radial  flow  component.  The  radial 
flow  component  near  the  surface  of  the  annulus  Is  of  particular  importance 
since  the  dynamic  pressure  of  this  flow  acts  against  the  blades  to  rotate 
the  annulus,  and  to  dislodge  and  expel 1  hydrometeors  from  the  annulus  sur¬ 
face. 

There  are  certain  critical  places  where  it  Is  important  to  know 
the  dynamic  pressure  on  the  surface  of  the  Instrument.  Especially  impor¬ 
tant  are  the  gaps  between  rotating  and  nonrotating  parts  at  which  the  interior 
of  the  instrument  is  exposed  to  the  environment.  With  reference  to  Figure 
5,  such  locations  are  the  hub  junction  with  the  annulus,  at  point  (x/r2  = 

0,  r/r2  *  2/3),  and  the  gap  between  the  annulus  and  the  downstream  housing, 
at  point  (x/r2  *  e,  r/r?  *  1).  At  the  hub  junction  the  x  component  of  flow 
is  zero,  and  extrapolation  of  calculated  velocities  shows  that  the  radial 
component  also  is  zero.  Thus,  we  have  stagnation,  and  the  dynamic  pressure 
relative  to  ambient,  aP$,  at  this  point  Is 

Aps  s  ?  p  v-  “  45’5 

where  p  Is  air  density  and  *  100  m  s"1  Is  free  stream  air  speed  (l.e., 
the  speed  of  the  airplane). 

The  closest  on-body  point  to  the  junction  of  the  annulus  with 
the  housing  for  which  we  have  a  velocity  calculation  Is  (x/r2  *  0.2, 
r/r2  *  1.0).  Here  we  have  vx/V„  *  2.4,  vr/V„  ■  0.  The  pressure 
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Figure  6.  Digital  description  of  the  disc  collector  with 
parabolic  hub. 
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Figure  7a.  Trajectories  of  water  drops  of  radius  0.5  ym  to 
the  disc  with  parabolic  hub.  These  trajectories  are  essentially 
coincident  with  flow  stream  lines.  Free  stream  air  speed  is 
100  *  s  .  r2  “  ^  CB1, 
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Figure  7b.  Trajectories  of  water  drops  of  radius  5  ym  to  the  di 
with  parabolic  hub.  E  is  impaction  efficiency. 
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Figure  7d.  Trajectories  of  water  drops  of  radius  35  to  the 
disc  with  parabolic  hub. 
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Figure  8.  Flow  field  near  the  annulus  surface  of  the  disc  with 
parabolic  hub.  (See  Figs.  5  and  6.) 
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relative  to  ambient  lb 


AP  =  lp  vf  [l  -  (vx/Vj2]  =  -63.6  mb. 

Thus  at  the  hub  junction  we  have  dynamic  pressure  acting  to  force  water 
into  the  interior  of  the  instrument,  whereas  at  the  housing  junction  we 
have  a  suction  outward  from  the  interior.  These  facts  have  been  accounted 
for  in  the  design  of  the  instrument  housing. 

DISC  WITH  PARABOLIC  HUB  AND  TRIANGULAR  RING 

A  further  modification  to  the  WCM  was  made  by  addition  of  a  ring 
with  triangular  crossection  upstream  of  the  outer  edge  of  the  rotating 
annulus  as  shown  in  Figure  10.  The  purpose  of  the  ring  is  to  supply  heat 
to  melt  ice  hydrometeors  and  to  shade  the  outer  edge  of  the  annulus  from 
direct  impaction.  As  will  be  seen  below,  the  ring  also  serves  to  increase, 
and  make  more  nearly  constant,  the  radial  airflow  component  near  the  annulus 
surface. 

We  also  consider  the  effect  on  the  flow  of  the  presence  of  the 

torque  blades.  As  indicated  in  Figure  11a,  the  blades  are  positioned  at 

12°  intervals  around  the  annulus  (there  are  30  blades),  and  the  width  of 

each  blade  covers  a  5°  arc  of  the  annulus  circle.  The  center  line  from  end- 

to-end  of  each  blade  is  inclined  at  an  angle  of  41.5°  to  the  radial  direc- 

* 

tion,  the  anqle  being  measured  at  the  outer  edge  of  the  annulus. 

Figure  lib  shows  the  projection  of  the  blades  in  the  crossection 
plane.  In  the  radial  range  r  *  2.0  to  2.5  cm  the  blade  height  is  1  mm; 
it  then  tapers  upward  to  a  height  of  5  mm,  which  height  is  reached  at 
r  =  2.7  cm,  and  then  remains  constant  to  r  =  r^  =  3  cm.  Since  the  annulus 


*This  corresponds  to  angle  e  of  Fig.  21. 
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a.  Plan  view  of  the  stationary  annulus  showing  the 
torque  blades. 


as  .•  •»  ••  •» 

r/rj 


b.  Crossectton  profile  of  the  disc  collector  with 
parabolic  hub  and  triangular  ring  showing  a 
projection  of  the  torque  blades  on  the  plane 
of  the  drawing. 


c.  Plan  view  of  the  annulus  showing  the  torque  blades 
approximately  as  seen  by  the  radial  flow  when  the 
annulus  Is  rotating  at  equilibrium  speed. 


Figure  11.  Torque  blade  geometry  for  the  water 
content  meter.  Drawings  are  to  proportional  scale 


ii 


and  blades  rotate  at  very  close  to  equilibrium  speed,  the  air  that 
passes  through  the  blades  sees  a  configuration  approximately  like  that 
of  Figure  11c.  We  have  included  the  tall  portion  of  the  blades,  as 
outlined  by  the  broken  lines  in  Figures  lib  and  11c,  in  the  digital  de¬ 
scription  (Fig.  12). 

Figures  13  and  14  show  the  flow  field  and  radial  component  of 
flow  near  the  surface  of  the  annulus.  The  plane  of  the  calculation  is 
midway  between  two  torque  blades.  The  presence  of  the  ring  and  blades 
is  seen  to  increase  the  radial  flow  speed  and  make  It  more  nearly  constant 
in  the  x  axis  direction  (cf.  Fig.  9). 

Figure  15  shows  the  variation  of  radial  flow  component  with 
r/r2  close  to  the  surface  of  the  annulus  (x/r2  =  -0.02).  Note 
that  the  curve  is  approximately  linear  out  to  r/r2  =  0.85. 

(A  similar  calculation  for  a  smaller  triangular  ring  and  slightly 
different  blade  geometry  also  was  calculated,  but  that  design  was  supplanted 
by  the  one  described  above,  and  its  description  would  add  nothing  of  value 
here.  ) 

HYDROMETEOR  SPECTROMETER 

At  the  time  of  this  writing  the  hydrometeor  spectrometer  (HS)  is 
in  a  preliminary  stage  of  development,  and  the  geometry  discussed  here  is 
tentative.  Similar  to  the  WCM,  the  measurement  surface  of  the  HS  is  a 
rotating  annulus  that  is  oriented  normal  to  the  free  stream  flow.  However, 
instead  of  having  a  stationary  hub,  the  central  portion  of  the  HS  is  open 
such  that  flow  is  allowed  through  it.  A  thin  wall,  called  the  "flow  splitter", 
rises  normal  to  the  annulus  surface  around  its  Inner  edge.  Thus,  we  have 
stagnation  at  the  inner  edge  of  the  annulus. 
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on  off  th« 
ades. 


Figure  12.  Digital  descrlptl 
triangular  ring  and  torque  bl 


Figure  13.  Flow  field  near  the  annulus  surface  of  the  WCM  with 
parabolic  hub,  triangular  ring  and  torque  blades.  The  plane  of 
the  calculation  Is  midway  between  two  torque  blades. 


25 


Figure  14.  Contours  of  constant  radial  flow  component,  v^V^,  near 

the  annulus  surface  of  the  WCM  with  parabolic  hub,  triangular  ring 
and  torque  blades.  The  calculation  plane  Is  midway  between  two  torque 
blades. 
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Figure  15.  Radial  flow  component  at  distance  x/r2  =  -0.02  from  the 
surface  of  the  WCM  annulus. 
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Inner  and  outer  radii  of  the  annulus  are,  respectively,  r^  =  5  cm 
and  r2  -  7  cm.  To  obtain  a  crude  approximation  to  the  flow  field  near  the 
surface  of  the  annulus,  the  geometry  shown  in  Figure  16  was  used  to  develop 
a  digital  description  of  the  annulus.  The  flow  splitter  is  not  actually  in¬ 
cluded  in  the  geometry;  instead,  the  inner  edge  of  the  annulus  is  extended 
inward  from  r/r2  =  0.714  to  0.5.  The  value  0.5  was  chosen  because  we 
estimated  that  the  stagnation  point  on  an  annulus  surface  of  these  dimensions 
should  be  near  the  desired  value  of  r/r2  =  0.714.  Figure  17  shows  the  digi¬ 
tal  description  used. 

Figures  18  and  19  show  the  flow  field  near  the  annulus  surface 
and  contours  of  the  radial  component  of  flow.  Stagnation  occurs  at  r/r?  * 
0.64  instead  of  0.714,  and  the  zero  contour  bends  toward  the  symmetry  axis 
instead  of  rising  vertically  from  the  annulus  surface  as  it  should  for  an 
accurate  modeling  of  the  flow  splitter  surface.  Nevertheless,  these  figures 
provide  a  qualitative  picture  of  the  expected  flow. 

Figure  20  shows  the  variation  of  the  radial  component  of  flow 
very  close  to  the  annulus  surface  (x/r?  =  -0.02)  with  r/r2<  The  curve 
is  approximately  linear  out  to  r/r?  =  0.9. 
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Figure  18.  Flow  field  near  the  annulus  surface  of  the  annular 
cylinder  shown  In  Figs.  16  and  17. 
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Figure  20.  Radial  flow  component  at  distance  x/r  =  -0.02  from  the  annulus  surface 
of  the  annular  cylinder  shown  in  Figs.  16  and  17. 


CALCULATION  OF  TORQUE  ON  THE  TORQUE  BLADES 


OUTLINE  OF  THE  ANALYSIS 

The  annuli  of  both  the  water  content  meter  (WCM)  and  hydrometeor 
spectrometer  (HS)  are  driven  rotationally  by  vanes,  called  torque  blades, 
that  are  mounted  at  an  angle  to  the  radial  flow  across  the  annular  surfaces 
as  shown  in  Figure  11  for  the  WCM.  Here  we  simplify  the  geometry  to  an 
annulus  with  a  single  blade.  As  shown  in  Figure  21,  the  blade  is  mounted 
at  an  angle  to  the  flow  such  as  to  cause  clockwise  rotation  of  the  annulus 
at  angular  speed  u. 

In  our  analytical  treatment  we  make  the  following  restrictions 
and  approximations: 

1.  The  blade  is  thin,  straight  and  of  constant  height,  h. 

2.  Radial  air  speed  across  the  annulus,  vr,  is  constant  with 
distance  from  the  surface  in  the  negative  x  axis  direction. 

3.  vf.  is  either  constant  with  r,  or  increases  linearly  with  r 
beginning  at  vr  *  0  at  r  =  r^ .  (See  Fig.  21  for  notation.) 

4.  The  annulus  is  rotating  at  or  near  to  its  equilibrium  speed, 
and  the  blade  offers  negligible  obstruction  to  the  flow. 

While  none  of  these  assumptions  are  in  accord  with  actuality,  the  last  is 
probably  the  most  unsatisfactory  from  a  theoretical  point  of  view.  This 
is  because,  as  indicated  in  Figure  21,  the  angle  a  between  the  radial 
direction  and  the  blade  varies  across  the  length  of  the  blade,  and,  as 
we  shall  see  below,  even  at  equilibrium  there  is  a  distribution  of  nonzero 
dynamic  pressures  across  the  blade.  Consequently,  there  will  be  a  reaction 
on  the  air  which  must  distort  its  flow. 

The  remainder  of  the  restrictions  can  be  relaxed,  and  more  exact 
calculations  can  be  done,  if  numerical  rather  than  analytical  integration 
methods  are  used. 
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DIFFERENTIAL  TORQUE 


We  begin  by  computing  the  dynamic  pressure  exerted  on  the  blade 
at  a  radial  distance  r  from  the  center  of  rotation.  At  this  point,  the 
radial  axis  makes  an  angle  a  with  the  blade.  The  blade  (and  annulus)  are 
rotating  at  speed  u. 

From  Figure  22,  we  see  that  in  a  frame  of  reference  on  the  ro¬ 
tating  blade,  air  flow  vector  is  v,  and  the  component  of  v  normal  to  the 
surface  of  the  blade  is  vp.  Therefore,  the  dynamic  pressure  on  the  blade, 
aP,  is 


aP  *  h  p  vp2 

It  is  easy  to  show  from  Figure  22  that 
Vp  *  s(-T  cos  a  +  J  sin  a) 


(2) 


(3) 


where 


s  *  vr  sin  a  -  wr  cos  a  (4) 

and  t  and  J  are  unit  vectors  in  the  directions  of  the  -u  x  r  and  r  axes. 

Consider  an  element  of  blade  area  of  differential  length  dL  and 
height  h  centered  at  our  point  of  interest.  The  differential  pressure  force 
on  this  element  of  area  Is 

d?  »  *sP  h|vp|vp  dL  .  (5) 

(The  use  of  | v p|  vp  instead  of  vp2  ensures  that  d t  has  the  proper  direction.) 
The  differential  torque,  dt,  on  the  blade  is 
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Figure  22.  Vector  diagram  that  shows  the  relationship  of  the  radial 
flow  vector,  vr,  the  flow  vector  from  the  perspective  of  the  rotating 
annulus,  v,  and  the  flow  component  normal  to  the  blade  surface,  vp. 
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d Y  =  rXd? 


(6) 


where  r  *  }r»  and  substitution  of  eqs.  (3)  and  (5)  into  (6)  qives 
dt  *  t  H  phr  cos  a|s|  s  dL 

where  £  is  the  unit  vector  normal  to  the  plane  of  Figure  22.  It  is  con¬ 
venient  at  this  point  to  drop  the  vector  notation  and  express  r  and  dL  in 
terms  of  a.  To  this  purpose  we  use  the  law  of  sines  to  derive  the  relation 


(see  Fig.  21) 

r  =  r2  sin  8  esc  a  (7) 

from  which  we  iitmediately  obtain 

rjsin  (<f>  +  s)  *  r2s1n  6  .  (8) 

Also,  it  is  easy  to  prove  that 

2 

dL  =  -r2sin  8  CSC  ada  ,  (9) 

and  from  this  we  obtain 

dT  *  Phr|s1n2e  |s|  s  da  (10) 

sin  a 


which  Is  our  basic  differential  torque  equation  in  dimensional  form. 

To  express  the  differential  torque  equation  in  nondimensional 
form,  we  define  a  nondimensional  rotation  speed, 

n  ■  «r2/vr  (11) 

where  *>  Is  rotation  speed  with  dimensions  radians-time 
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a  nondimensional  s(a) 


S(a)  =  sina  -  asing  etna  (12) 

and  a  nondimensional  torque, 

T  *  T/(>*h  vr2r22)  (13) 

where  T  is  torque  with  dimensions  force-length, 
in  terms  of  which  eq.  (10)  becomes 

cfT  =  -sin2e|S(a)|  S(a)  cosa  csc\  da,  (14) 

which  is  our  basic,  nondimensional ized  differential  torque  equation. 

It  is  apparent  that  the  sign  of  the  differential  torque  depends 
on  the  sign  of  S(a).  Specifically,  if  S(a)  >  0  the  torque  boosts  rotation, 
whereas  if  S(a)  <  0  the  converse  is  true.  At  equilibrium  or  near  equilibrium 
rotation  speed,  there  will  be  a  point  along  the  length  of  the  blade  at  which 
S(a)  is  zero;  let  this  point  be  characterized  by  the  angle  aQ. 

To  determine  total  torque  on  the  blade,  eq.  (14)  must  be  integrated 
from  a  =  e  +  ♦  to  e  as  implied  by  Figure  21.  For  actual  flows  and  blade 
geometries,  this  would  have  to  be  done  numerically.  Nevertheless,  it  is 
worthwhile  to  study  two  simple  cases  for  which  the  integration  can  be  done 
analytically.  First  we  consider  the  case  with  vp  constant,  and  then  we 
consider  linear  variation  of  vp  with  r.  In  both  cases  we  take  h  to  be  con¬ 
stant. 

TOTAL  TORQUE  WITH  CONSTANT  RADIAL  FLOW 

With  both  vp  and  h  constant  eq.  (14)  can  be  expressed  in  the  form 
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05) 


2  2 
dT  =  sin  e  etna  (1  -  nsine  etna  csca)  da 


and  we  obtain  the  Indefinite  integral 
Y(a,fl,g)  =  J |d  T| 


=  sin20 


m(sina)  +  nsine  £ctna  csca  +  tn(tan£)J 


2  2 

-  (fc)sin$  ctn  a) 


(16) 


As  noted  above,  dT  changes  sign  at  a  =  ao.  Therefore,  we  must  determine 
If  aQ  lies  in  the  range  $  +  $  to  $,  and  if  so,  the  integration  of  equation 
(15)  must  be  done  in  two  parts:  from  a  =  <t>  +  8  to  aQ,  and  from  aQ  to  8, 
with  the  appropriate  sign  being  taken  for  each  part.  Accordingly,  we  have 
the  following  situations: 

'•  <  0,  we  calculate  <*o  from  the  equation 


aQ  =  tan”  ‘V  +  (2/nsin8)2  -  l) 


(17) 


and  calculate  the  total  torque  from 


T  *  Y(8,n.e)  +  y(*  +  B.n,e)  -  2Y(ao,n,e) 


(18) 


2. 

we  calculate  the  total  torque  from 


T  *  y(4  +  8,«,e)  -  Y(B.n.B) 


(19) 


but  change  the  sign  of  the  result  if  S(e)  <  0. 
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For  given  geometry  there  is  an  equilibrium  rotation  speed, 
fte,  at  which  there  is  zero  torque  on  the  blade.  We  have  calculated  ng 
as  a  function  of  e  for  the  WCM  geometry  (i.e.,  r-|  =  2  cm,  r2  =  3  cm) 
by  Newton-Raphson  iteration  as  explained  in  Appendix  B,  and  the  result 
is  plotted  in  Figure  23.  Notice  that  8  is  limited  to  the  range 

0  <  6  <  8  (20) 

where  from  eq.  (7)  and  Figure  21  we  see  that 

6max=si"'1(7Js1n90“)  (2,) 

For  the  WCM  geometry  we  have  8  =  41.81°,  and  from  Figure  23  we  have 

a  =  1.6915  which  corresponds  to  uk  =  4511  rad  s"1  (4.307  X  104RPM)  for 
e  i  e 

vr  =  80  m  s"  .  (See  Fig.  15.) 

Figure  24  shows  T  vs  a  for  B  *  33.094°  U  +  8  *  54.987°)  and 
ag  =  1.044.  Taking  vr  =  80  m  s"1,  this  latter  corresponds  to  u>e  =  2784 
rad  s'1  (2.658  X  104  RPM). 

Also  of  interest  is  the  distribution  of  dynamic  flow  pressure 
across  the  length  of  the  blade.  For  this  we  define  a  dimensionless  dynamic 
pressure  (see  eq.  (2),  ff. ) 

AP/0spVr2)  =  |S(a)|  S(a)  (22) 

where  S(a)  is  given  by  eq.  (12).  This  is  plotted  for  the  WCM  geometry 

with  B-e  in  Figure  25  as  a  function  of  o  for  several  values  of  a. 
max 

Figure  21  shows  the  geometric  relationship  between  angle  o  and  distance 
along  the  blade.  Quantitatively  this  relationship  is 
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Figure  23.  Nond linen si one 1  equilibrium  rotation  speed  vs  blade  angle 
8  for  the  WCM  with  constant  radial  flow. 
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33.09 


Figure  24.  Nondimensional  torque  (eq.  (13))  vs  nondimensional  rotation  speed  (eq.  (11)) 
for  the  WCM  with  blade  angle  e  *  33.09°.  Radial  flow  is  constant. 
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Figure  25.  Dynamic  pressure  distributions  across  the  length  of 
a  MCM  blade  oriented  at  maximum  angle,  for  several  rotation 

speeds.  Radial  flow  Is  constant.  The  heavy  curve  Is  for  equilibrium 
rotation  speed.  Positive  pressure  boosts  rotation  speed;  negative 
pressure  retards  rotation. 
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where  t  is  the  distance  along  the  blade  from  its  inner  end,  and 

L  =  r2cos0  -  rjeosU  +  e).  (24) 

For  the  WCM  with  8  =  ft.  we  have  l  =  3  cosa  cm  and  L  =  2.236  cm. 

ma  X 

At  the  equilibrium  rotation  speed,  the  crossover  from  boost  to  drag  occurs 
at  aQ  =  54.26°  (t/L  =  0.784).  Thus,  it  is  apparent  that  torque  will  be 
increased  (reduced)  by,  for  example,  decreasing  (increasing)  blade  height 
toward  the  outer  edge  of  the  blade  relative  to  the  inner  edge.  Of  course 
such  a  change  also  will  affect  equilibrium  rotation  speed  in  a  complex 
manner  that  would  have  to  be  determined  by  further  calculation  or  experiment. 
It  should  also  be  borne  in  mind  that  the  assumption  of  constant  vr  is  un¬ 
realistic  (cf.  Fig.  15). 

Next,  we  consider  the  more  realistic  assumption  of  linear  increase 
of  vr  with  r,  but  use  the  hydrometeor  spectrometer  for  our  case  study  instead 
of  the  water  content  meter. 

TOTAL  TORQUE  WITH  LINEAR  VARIATION  OF  RADIAL  FLOW 

In  this  case  we  take  vf  to  be  zero  at  the  inner  edge  of  the  annulus 
(i.e. ,  at  r  =  r^  and  vf  to  be  2vr  at  the  outer  edge  of  the  annulus  (i.e., 
at  r  =  r2),  where  vr  is  the  mean  radial  speed  across  the  annulus,  and  we 
assume  linear  variation  of  vp  between  these  extremes.  Thus  we  have 


45 


Using  eqs.  (7)  and  (8)  this  becomes 


2  v 


vr  =  -jp  sing  £csca  -  csc($  +  g)J  , 

and  from  this  result  plus  eqs.  (4),  (11)  and  (12)  we  derive 
S(a)  =  —  =  |  [sing  -  (1  -  R)sinaJ  -  nsing  etna 


where 


-  r2  ~  rl  _  sinU  +  g)  -  sing 
r0  sin(<|»  +  g) 


Corresponding  to  eq.  (16),  we  have 
D{a,a,B.R)  (list )2ct„2a  nRctrtt[«.J.  nR  etna)  -2 

+  2 

As  is  the  case  with  constant  vf,  we  have  two  situations: 

1.  If  <  0,  the  dynamic  pressure  reverses  sign  at  angl 

a  ,  where  g  <  a  c  $  +  g,  and  we  compute  total  torque  from 


)2 


*n(sin  a) 


(V)(^) 


4  csca  -  nR  etna  csca  +  tn(tan 


t>] 


(30) 


T  =  D(8,fl,e,R)  +  D( 4>  +  B.n.B.R)  -  2D(ao>£2,B,R) 
but  reverse  the  sign  of  the  result  if  2sine  >_  a  cosp. 

2.  If  >  0,  we  compute  total  torque  from 

T  =  DU  +  B.n.B.R)  -  D(B,Q,B,R)  (31) 

but  reverse  the  sign  of  the  result  if  2sine  £  n  cose. 

In  the  normalization  of  the  torque,  eq.  (13),  we  replace  vr  by  vr  for  this 
case. 

Calculation  of  aQ  is  not  as  straight  forward  as  for  the  case  with 
vr  constant.  If  we  set  eq.  (27)  to  zero  we  obtain 

n  sine  etna  +  2(- I- R)sina  =  2^^- 

o  K  o  R 

which  can  be  expressed  as 

Vl  +  C2  Me  -  1)  +  (1  -  R)/sine  =  0  (32) 

where  e  =  ctnaQ.  This  is  a  quartic  in  e,  with  four  roots,  of  which  at  least 
three  may  be  real.  Equation  (32)  is  easiest  solved  by  Newton-Raphson  itera¬ 
tion,  with  care  being  taken  to  obtain  the  relevant  root. 

The  hydrometeor  spectrometer  (HS)  is  used  as  the  particular  case 
to  study  these  equations.  Preliminary  specifications  for  the  HS  are 
r,  =  5  cm,  r,  =  7  cm  for  which  we  obtain  (eq.  (21))  e_.v  =  45.585°. 

I  L  mOX 

Figure  26  displays  T  vs.  n  for  &  ~  s  v  and  &  =  40°.  From  the  figure 

nlaX  j 

we  find  £l_(B_.v)  =  1.7825.  If  we  take  v„  =  50  m  s“  (see  Fig.  20),  then 

c  UldX  j  I 

a>e  =  1273  rad  s"  (1.216  X  10  RPM).  Also  we  note  that  torque  is  comparatively 

insensitive  to  rotation  speed  in  the  neighborhood  of  equilibrium  rotation. 

Pressure  distributions  across  the  length  of  the  blade  are  computed 

from  eq.  (22)  with  S(a)  given  by  eq.  (27).  Results  are  shown  in  Figure  27 

for  6  =  B  and  6  *  40°.  Notice  for  the  equilibrium  rotation  speeds  that 
max 
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26.  Nondlmenslonal  torque  \eq.  (13)}  vs  nondlmenslonal  rotation  spee< 
1))  for  the  HS  annulus  with  two  values  of  blade  angle  6.  Radial  flow 
linearly  with  r. 


a  (4tgr««s) 


Figure  27a.  Dynamic  pressure  distributions  across  the  length  of  an 
HS  blade  oriented  at  maximum  angle  3max  for  several  rotation  speeds. 

Radial  flow  varies  linearly  with  r.  The  heavy  curve  Is  for  equilibrium 
rotation  speed.  Positive  pressure  boosts  rotation  speed;  negative 
pressure  retards  rotation. 
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Figure  27b.  Dynamic  pressure  distributions  across  the  length  of  an 
HS  blade  oriented  at  maximum  angle  &  ■  40°  for  several  rotation  speeds. 

Radial  flow  varies  linearly  with  r.  The  heavy  curve  Is  for  equilibrium 
rotation  speed.  Positive  pressure  boosts  rotation  speed;  negative  pressure 
retards  rotation. 
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positive  dynamic  pressure,  which  boosts  rotation,  occurs  toward  the  outer 
edge  of  the  annulus,  while  there  is  drag  toward  the  inner  edge.  This  is 
precisely  opposite  to  the  situation  found  for  the  WCM  meter,  for  which  v 

r 

is  taken  to  be  constant  (cf.  Fig.  25).  Additional  comparisons  are  made 
in  the  next  section. 

SUMMARY  AND  DISCUSSION 

We  have  derived  an  approximate  method  for  calculation  of  dynamic 
pressure  distributions  across  the  length  of  water  content  meter  (WCM)  and 
hydrometeor  spectrometer  (HS)  torque  blades,  and  total  torque  on  the  blades. 
Calculation  of  total  torque  requires  solution  of  a  differential  equation, 
which  has  been  done  analytically  for  two  simple  cases: 

1.  Constant  radial  flow  across  the  annulus  face. 

2.  Linear  variation  of  radial  flow  across  the  annulus  face. 

Blade  height  is  taken  to  be  constant  for  both  cases.  The  first  case  is 
applied  to  the  WCM  geometry,  and  the  second  to  the  HS  geometry. 

Graphs  of  dimensionless  total  torque  vs  dimensionless  rotation 
speed  for  the  two  cases  are  similar  (cf.  Figs.  24  and  26).  On  the  other 
hand,  comparison  of  the  dynamic  pressure  distributions  for  equilibrium 
rotation  for  the  two  cases  (Figs.  25  and  27),  shows  reverse  distributions 
of  pressure  signs  across  the  lengths  of  the  blades  for  the  two  cases. 

While  the  geometries  are  somewhat  different  (i.e.,  »  o.333  while 

Rh$  *  0.286  (see  eq.  (28)),  the  major  difference  is  in  the  radial  flow  as 
stated  above,  and  this  shows  a  potential  for  high  sensitivity  to  change  in 
radial  flow. 

We  have  not  Investigated  the  effect  of  variation  of  blade  height 
on  total  torque  or  on  equilibrium  rotation  speed.  However,  the  pressure 
distribution  curves  Indicate  that  at  least  the  total  torque  would  be 
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sensitive  to  such  variation.  In  particular  Figure  27,  which  is  for  the 
more  realistic  case,  indicates  that  total  torque  would  Increase  with  in¬ 
crease  of  blade  height  toward  the  outer  edge  of  the  annulus.  The  effect 
of  this  change  on  equilibrium  rotation  speed  Is  not  Immediately  obvious. 

In  addition,  we  have  not  studied  the  effect  of  variation  of  flow 
in  the  direction  normal  to  the  annulus  face.  Figures  14  and  19  show  that 
this  variation  in  flow  can  be  significant  even  without  consideration  of 
the  boundary  layer.  Though  an  accurate  study  would  require  calculation  of 
the  boundary  layer,  useful  information  could  be  obtained  by  extension  of 
the  methods  used  in  this  study. 

Calculations  could  be  done  with  realistic  variation  of  both  blade 
height  and  flow  by  numerical  solution  of  eq.  (10).  Sensitivities  to  these 
as  well  as  other  parameters  could  be  studied,  and  the  results  then  could  be 
used  to  determine  the  torque  blade  structure  that  would  come  closest  to  pro¬ 
viding  the  desired  balance  of  rotation  speed  against  torque  sensitivity. 


HYDROMETEOR  IMPACT 


INTRODUCTION 


Here  we  are  concerned  with  effects  of  hydrometeor  impact  other 
than  those  which  are  used  for  measurement,  but  which  may  affect  instrument 
performance.  We  consider  lateral  jetting  of  impacting  water  drops,  splash 
and  bounce,  erosion,  and  stress  limits  on  torque  blades  imposed  by  collisions 
with  water  drops.  We  make  use  of  the  literature  and  some  theoretical  cal¬ 
culations;  no  new  data  are  reported  here. 

HIGH  SPEED  IMPACT  OF  WATER  DROPS:  LATERAL  JETTING 

As  a  large  water  drop  impacts  at  high  speed  on  a  rigid  surface, 
a  shock  wave  is  reflected  back  into  the  drop  as  shown  in  Figure  28.  Com¬ 
pressed  water  is  contained  between  the  shock  front  and  the  rigid  surface 
until  the  speed  of  the  shock  front  exceeds  the  speed  of  growth  of  the  inter¬ 
face,  after  which  the  shock  front  detaches  from  the  surface.  Immediately 
after  detachment,  the  pressure  is  relieved  by  a  jetting  of  water  at  high 
speed  along  the  rigid  surface.  The  problem  at  hand  Is  to  estimate  whether 
or  not  a  significant  amount  of  water  might  be  jetted. 

According  to  Heymann^  the  critical  contact  angle,  4>c,  at  which 
jetting  begins  for  an  Impact  speed*  of  100  m  s’  Is  about  2.5°.  However, 
Brunton  and  Camus ^  assert  that  Heymann's  *c  values  are  much  lower  than 


3.  F.  J.  Heymann,  "High-Speed  Impact  Between  a  Liquid  Drop  and  a  Solid 
Surface,"  J.  Appl.  Phys.  40,  5113  (1969). 

4.  0.  H.  Brunton  and  J.  J.  Camus,  "The  Flow  of  a  Liquid  Drop  During 
Impact,"  Proceedings  of  the  3rd  International  Conference  on  Rain 
Erosion  and  Associated  Phenomena.  A.  A.  Fyall  and  R.  B.  King,  edi¬ 
tors,  (Elvetham  Hall,  11-13  Aug.  1970).  pp.  327  ff. 

*A1rplane  cruising  speed  Is  taken  here  to  be  about  100  m  s’1. 


53 


Figure  28.  Development  of  a  reflected  shock  front  end  compressed 
Miter  (hitched  irei)  during  Impact  on  a  rigid  surface  of  a  large 
Miter  drop. 


observed,  and  they  report  values  between  12°  to  20°.  (They  also  report 
jetting  speeds  of  slightly  higher  than  1000  m  s"  for  100  m  s“  impact 
speed. ) 

Since  jetting  occurs  only  long  enough  to  relieve  the  compression 
behind  the  shock  front  we  can  assume  that  a  volume  of  water  less  than  that 
trapped  behind  the  shock  front  at  <j>  =  <j>c  will  be  jetted.  If  we  assume  no 
distortion  of  drop  or  shock  front,  the  volume  of  water  between  the  shock 
front  and  the  rigid  surface  is  simply  a  spherical  segment  of  one  base  (see 
Fig.  28),  and  the  fraction  of  the  total  drop  volume  involved,  f,  is 

f  3  h  (1  -  cos<f>c)2(2  +  cos<|>c).  (33) 

For  the  several  values  of  $c  mentioned  above  we  have: 


♦c 

f 

2.5° 

6.8 

X 

o 

1 

12.0° 

3.6 

1 

o 

X 

20.0° 

2.7 

1 

o 

X 

These  numbers  indicated  that  jetting  is  not  likely  to  be  a  serious  problem, 
though  we  must  caution  that  since  we  have  Ignored  the  effects  of  rotation, 
the  applicability  of  these  results  to  the  WCM  and  HS  is  questionable. 

LOW  SPEED  IMPACT  OF  WATER  DROPS 

Normal  Incidence 

Mutchler  and  Hansen ^  present  rather  comprehensive  experimental 
data  on  impact  on  stationary  surfaces  of  large  water  drops  at  terminal 

5.  C.  K.  Mutchler  and  L.  M.  Hansen,  "Splash  of  a  Water  Drop  at  Terminal 
Velocity,"  Science  169,  1311  (1970). 
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settling  speeds.  Their  most  relevant  information  is  summarized  as  follows: 

a.  Impact  on  dry,  smooth,  hard  surfaces  produces  no  splash;  the 
drop  merely  spreads  horizontally. 

b.  Maximum  splash  effects  occur  for  surface  water  depths  of  about 
1/4  to  1/3  of  the  drop  diameter. 

c.  Maximum  splash  height  is  about  6  times  the  drop  diameter. 

Using  the  data  of  Gunn  and  Kinser^,  we  estimate  that  maximum  impact 
speed  was  about  9  m  s“  . 

Oblique  Incidence 

The  only  data  found  are  in  two  papers  by  Taylor^7,8^  on  impact 
of  identical  water  jets  at  various  angles  of  incidence.  Sketches  of 
Taylor's  results  are  shown  in  Figure  29. 

After  Impact,  the  jets  merge  to  form  the  thin  sheets  shown  in 
the  figure,  where  the  sheets  lie  In  the  plane  which  bisects  the  angle  be¬ 
tween  the  jets.  (This  plane  can  be  replaced  by  a  rigid  surface  onto  which 
one  of  the  jets  is  Impacting  at  an  angle  of  0.)  Since  friction  and  surface 
effects  are  relatively  small  In  the  region  of  Impact  of  the  jets,  momentum 
conservation  requires  that  flow  speeds  In  the  sheets  be  essentially  equal 
to  the  jet  speed.  Impact  pressure  effects  cause  the  sheet  thickness  to 
vary  with  angle  $.  Then,  when  surface  tension  pressure  at  the  edge  of 
the  sheet  balances  dynamic  flow  pressure  (i.e. ,  when  the  Weber  number  at 
the  sheet  edge  reaches  unity),  the  sheet  br.-aks  up  into  tiny  droplets 


6.  R.  Gunn  and  G.  D.  Kinzer,  "The  Terminal  Velocity  of  Fall  for  Water 

Droplets  In  Stagnant  Air,"  0.  Meteor.  6,  243  (1949). 

7.  G.  I.  Taylor,  "Formation  of  Thin  Flat  Sheets  of  Water,"  Proc.  Roy. 
Soc.  (London)  A259,  1  (1960). 

8.  G.  I.  Taylor,  "Oblique  Impact  of  a  Jet  on  a  Plane  Surface,"  Phil. 
Trans.  Roy.  Soc.  (London)  A260,  96  (1966). 
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which  fly  off  tangentially  to  the  edge  of  the  sheet.  This  limits  the 
breadth  and  length  of  the  sheets  as  indicated  in  the  figure. 

Jet  speeds  of  3.6  m  s-1(e  =  60°  and  45°)  and  52  ms"  (e  =  30°) 

were  used. 

MAXIMUM  REBOUND  DISTANCE  AND  RETURN  TIME 

To  calculate  maximum  particle  rebound  after  impaction  and  maxi¬ 
mum  return  time  to  the  surface  after  initial  impact  and  rebound,  we 
assume  perfectly  elastic  impact.  Thus,  we  ignore  fragmentation  and  in¬ 
elastic  effects.  Also,  we  ignore  the  annulus  rotation  and  the  boundary 
layer. 

We  use  the  results  of  the  particle  trajectory  calculations  to 
the  WCM  with  parabolic  hub  {Figs.  7)  to  determine  particle  impact  speeds, 
which  are  plotted  in  Figure  30. 

Air  flow  components  normal  to  the  annulus  surface,  as  computed 
by  the  Hess-Smith  code  for  the  body  shown  in  Figure  6  at  a  distance  of 
r  =  2.4  cm  from  the  WCM  axis,  were  used  in  the  particle  trajectory  calcula¬ 
tions.  This  flow,  v  /V  ,  was  found  to  be  closely  fit  by  the  equation 

n  -  c[l  ±  Vl  +  (0.215/c)2  ]  (34) 

where 

n  -  *og10  (vx/Vj 
C  s  H  [*°910  (-x/r2)  +  0.518] 

and  the  +  sign  is  taken  if  c  Is  minus  while  the  -  sign  is  taken  if  c 
is  positive. 

Initial  particle  rebound  speeds  were  taken  from  Figure  30,  and 
were  directed  in  the  minus  x  direction.  Air  flow  was  calculated  from  eq.  (34). 
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Calculated  water  drop  impact  speeds,  v  ,  on  the  disc  annulus  with  parabolic 


Trajectories  were  calculated  for  spheres  by  use  of  the  methods  described 
in  reference  1 ,  except  that  gravity  was  ignored.  Atmospheric  conditions 
for  3  km  altitude  in  the  U.  S.  Standard  Atmosphere  were  used: 

T  =  268.66°  K 

P  =  70121  Pa 

p  =  0.90925  kg  m"3 

n  =  1.6938  X  10'5  kg  m'1  s"1 

Free  stream  air  speed,  V  ,  was  taken  to  be  100  m  s-1.  Particle  densities  of  917 

3  °°  ••  3 

kg  m"  (ice),  and  100  kg  m  (snow)  were  used.  Results  are  listed  in  Table  1. 

EROSION 

Water  Drops 

(41 

According  to  Brunton  and  Camus'  ,  surface  erosion  in  the  form 
of  "small,  shallow  depressions  or  dimples"  begins  to  appear  on  stainless 
steel  after  about  1000  water  drop  impacts  at  125  m  s"  .  From  this,  and 
other  data  presented  in  the  Proceedings  of  the  Third  International  Con¬ 
ference  on  Rain  Erosion  and  Associated  Phenomena it  appears  that  at  the 
comparatively  low  impact  speeds  of  Interest  here  there  would  not  be  an 
immediate  erosion  problem.  However,  after  extended  use,  the  instrument 
surfaces  probably  would  become  roughened  by  erosion,  particularly  if 
constructed  of  a  soft  metal  such  as  aluminum. 

Parker  makes  an  observation  of  interest  in  his  addendum  to  refer¬ 
ence  4.  He  found  that  surface  roughness  causes  a  marked  decrease  in  water 
drop  retention  on  a  rotating  surface.  Moreover,  he  shows  sketches  of 
disintegration  of  an  impacted  water  drop  into  small  droplets  after  low- 
speed  impact  on  a  rotating  surface,  the  disintegration  presumably  being 
caused  by  surface  roughness. 


60 


TABLE  1 


MAXIMUM  PARTICLE  REBOUND  DISTANCES  AND  RETURN  TIMES 


a.  Particle  density  917  kg  m 


Initial 

Max. 

Time  at 

Particle 

Rebound 

Rebound 

Max.  Rebound 

Return 

Final 

Diameter 

Velocity 
(m  s-1) 

Distance 

Distance 

Time 

Velocity 
(m  s"  ) 

(pm) 

(cm) 

(s) 

(s) 

10 

-  63.0 

-  .403 

1 .485X10  4 

4. 38X1  O'*4 

13.2 

20 

-  74.0 

-  .921 

2.700 

7.26  3 

24.1 

40 

-  88.4 

-  2.116 

5.400  3 

1.38X10 

33.7 

70 

-  93.7 

-  4.252 

1.05  X10" 

2.64 

38.9 

100 

-  95.6 

-  6.469 

1.62 

4.04 

40.4 

200 

-  97.5 

-  13.93 

3.54 

8.88  , 

41.1 

400 

-  98.5 

-  28.86 

7.44  , 

1.87X10 

41.2 

700 

-  98.9 

-  51.08 

1.32  X10 

3.34 

41.3 

1000 

-  99.2 

-  73.28 

1.89 

4.80 

41.4 

2000 

-  99.5 

-146.7 

3.78 

9.63  . 

1.927X10' 

41.6 

4000 

-100.0 

-293.7 

7.56  , 

1.32  XI  O' 

41.8 

7000 

-100.0 

-511.8 

3.36 

41.9 

b.  Particle  density  100  kg  m'3 

Initial 

Max. 

Time  at 

Particle 

Rebound 

Rebound 

Max.  Rebound 

Return 

Final 

Diameter 

Velocity 

Distance 

Distance 

Time 

Velocity 

(pm) 

(m  s'1) 

(cm) 

(s) 

(s) 

(m  s'1) 

10 

-  63.0 

-  .076 

4.05X10'5 

3.48X10"4 

0.11 

20 

-  74.0 

-  .206 

8.10  4 

3.09 

5.56 

40 

-  88.4 

-  .462 

1.32X10" 

4.30 

15.2 

70 

-  93.7 

-  .811 

2.01 

5.94 

23.2 

100 

-  95.6 

-  1.15 

2.76 

7.70  3 

28.0 

200 

-  97.5 

-  2.14 

5.16 

1.37X10' 

33.5 

400 

-  98.5 

-  3.97 

9.60  , 

2.54 

36.4 

700 

-  98.9 

-  6.58 

1.63X10' 

4.24 

37.9 

1000 

-  99.2 

-  9.15 

2.28 

5.93  „ 

38.6 

2000 

-  99.5 

-  17.52 

4.44 

1.14X10 

39.7 

4000 

-100.0 

-  34.08 

8.76  2 

2.23 

40.5 

7000 

-100.0 

-  58.5 

1.51X10' 

3.83 

40.9 
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Ice  Particles 


The  question  here  is  whether  significant  erosion  is  to  be  expected 
from  impact  of  large  single-crystals  of  ice  at  impact  speeds  of  about  100 
m  s"1.  We  have  not  succeeded  in  finding  specific  data  with  which  to  answer 
this  question.  There  are,  however,  good  general  discussions  and  some  inter¬ 
esting  data  in  the  book  by  Engle^.  Here  we  summarize  some  of  Engle's 
more  pertinent  discussions. 

Impact  surface  damage  is  caused  by  the  immediate  effects  of  the 
collision  (primary  erosion),  and  also  by  secondary  collisions  of  fragments 
(secondary  erosion).  Thus,  the  dent  or  pit  caused  by  the  primary  erosion  is 
usually  surrounded  by  the  many  smaller  scars  of  the  secondary  erosion. 

Primary  erosion  damage  on  ductile  surfaces  is  usually  a  maximum  at  angles 
of  incidence  between  20-30°,  but  secondary  erosion  reaches  a  maximum  at 
90°  incidence. 

For  normal  incidence,  erosion  is  independent  of  particle  size 
over  a  considerable  range  of  sizes. 

It  has  been  found  that  weaker,  fragile  particles  sometimes  pro¬ 
duce  greater  erosion  than  stronger  particles.  This  is  caused  by  fragmertati . 
where  the  secondary  particles  remove  lips  formed  by  primary  impact. 

Summary 

At  the  impact  speed  of  about  100  m  s-1  under  consideration  here, 
it  is  apparent  that  water  drops  will  not  cause  significant  surface  erosion 
except  perhaps  after  extended  exposure.  In  the  case  of  large  ice  particles, 
however,  we  cannot  be  assured  that  erosion  will  be  insignificant,  expecially 
if  a  comparatively  soft  metal  such  as  aluminum  is  used.  Moreover,  we  can 


9.  P.  A.  Engel,  Impact  Wear  of  Materials  (Elsevier,  1976).  Chapter  5. 
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expect  that  surface  roughness  caused  by  erosion  will  alter  the  retention 
characteristics  of  the  instrument  surfaces,  and  thereby  affect  instrument 
performance. 

IMPACT  STRESS  ON  THE  TORQUE  BLADES 

We  are  concerned  here  with  the  force  of  impact,  and  resulting 
impact  stress,  of  large  water  drops  on  the  sides  of  the  torque  blades. 
Impact  is  caused  by  rotation  of  the  blades  into  the  water  drops.  The 
questions  to  be  answered  are:  how  much  impact  stress  is  generated  in 
the  blades;  which  candidate  construction  metals  are  satisfactory  and  which 
are  not;  and  how  thick  must  the  blades  be  to  reduce  the  stress  to  tolerable 
levels? 


Force  of  Impact  of  a  Large  Water  Drop 

A  plan  view  of  the  overall  annulus  and  blade  geometry  is  shown 
in  Figure  21,  and  details  of  the  WCM  blade  geometry  are  shown  in  Figure  11. 
From  Figure  21  and  eq.  (11),  it  is  apparent  that  the  maximum  impact  speed, 
v^  is 

v..  =  a>r2cos0  =  vrncose  .  (35) 

According  to  Heymann^  the  peak  impact  pressure  of  a  large  drop 
is  approximately 


where  pw  Is  water  density  and  Cw  is  speed  of  sound  in  water.  However,  this 
seems  to  be  localized  at  the  periphery  of  the  Impact  area,  while  at  the 
center  of  the  impact  area,  the  pressure  is  closer  to 
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Pcenter  “  pw^wvi 


For  this  calculation  we  will  use  a  simple  average  of  these  equations  for 
the  impact  pressure,  P^, 

pi  ■  •  <36> 

Critical  contact  angle,  4  ,  may  be  defined  as  the  contact  angle  (Fig.  28) 
at  peak  pressure.  From  Figure  28  the  contact  area  at  Deak  oressure,  A^. 
is 

A.  =  it  (a  sin  $c)2  ,  (37) 

where  a  is  drop  radius.  Thus  the  impact  force,  Fi  =  P^A..,  is  given,  on 
combining  eqs.  (35)-(37),  by 

F..  =  1.4k  PwCwa2vrnsin2<frccose  .  (38) 

Maximum  Stress  on  the  Blade 

Standard  engineering  methods  used  for  calculation  of  maximum 
stress  in  a  beam  can  be  used  here^10^.  Geometry  and  notation  are  shown  in 
Figure  31.  The  Impact  force  causes  the  blade  to  bend,  which  sets  up  stress 
forces,  S,  in  the  blade  which  act  parallel  to  the  surface  of  the  blade  as 
shown  in  the  Figure.  On  the  convex  (upper)  surface  of  the  bent  blade,  the 
stress  is  caused  by  stretching,  while  on  the  concave  (lower)  surface,  the 
stress  is  caused  by  compression.  At  the  surfaces,  the  stresses  are  at 
their  maxima  and  are  designated  as  Sm.  For  our  blade  of  height  h,  width 
L  and  thickness  s,  the  maximum  stress  is 


10.  E.  Hausmann  and  E.  P.  Slack,  Physics  (D.  Van  Nostrand,  3rd  Edition, 
1948).  Secs.  107-108. 
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,  _  M6 
\  "  21 


where 


(39) 


M  *  Fjh  (40) 

is  the  bending  moment,  and  I  is  the  moment  of  Inertia  of  the  blade  crossec- 
tion  area,  given  by 

I  *  L63/12  (41) 


for  a  blade  of  rectangular  crossectlon. 

Combining  eqs.  (38)-(41)  we  have,  finally, 


m 


8.4*  PwCwa2vrncoses1n2*c 


(41) 


For  the  WCM  we  can  use  %  ■  aB#x  and  ne  -  1.6915  (Fig.  23).  From 
Figure  14,  we  take  vr  *  80m  s~  as  a  reasonably  conservative  typical  value. 
Brunton  and  Camus^  report  that  #c  Is  in  the  range  12-20°;  we  use  ac  *  15° 
here.  Other  quantities  are  pw  -  1000  kg  nf  and  Cw  *  1500  m  s’  . 

From  Figure  11,  we  see  that  the  height  of  the  tall  portion  of 
the  WCM  blade  Is  5  mm.  To  allow  for  full  Impact  of  large  water  drops  near 
the  end  of  the  blade,  we  take  h  *  4  mm.  The  width  of  the  tall  portion  of 
the  blade  Is  roughly  L  ■  6  mm.  Thus,  in  terms  of  water  drop  radius  and 
blade  thickness,  we  obtain 


1.783  X  108(a/«)2  (N  i'2) 

1.783  X  10*(a/«)2  (dynes  cm’2) 

2.586  X  lO^a/a}2  (ptl) 


Equation  (42)  Is  plotted  In  Figure  32. 


(42) 
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MAXIMUM  STftCSS 


Figure  32.  Maximum  tensil  stress  vs,  blade  thickness  for  impact 
of  water  drops  of  radius  a  according  to  eq.  (42). 
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Yield  Strength  and  Tensil  Strength  of  Metals 

Tensil  and  yield  strengths  of  metals  are  determined  by  applying 
stretching  forces  to  bars  of  metal  and  measuring  the  linear  stretch  as  a 
function  of  applied  force.  The  tensil  stress,  S,  is  defined  as 

S  =  Force/original  crossectional  area  of  bar. 

This  is  plotted  against  strain,  s,  which  is  defined  as 

s  =  Increase  in  length  of  bar/original  length  of  bar. 

A  typical  stress-strain  curve  is  shown  in  Figure  33.  The  elastic  ranqe, 
over  which  the  bar  returns  to  Its  original  length  when  the  load  is  re¬ 
leased,  is  the  range  of  applicability  of  Hook's  law.  Yield  strength  is 
defined  to  be  the  tensil  stress  at  a  point  at  which  the  permanent  set 
(i.e.,  permanent  stretch)  is  0.2%  of  the  original  bar  length.  Tensil 
strength  Is  the  maximum  tensil  stress  on  the  curve.  (Tensil  strength  is 
also  called  "ultimate  tensil  stress"  or  "ultimate  strength.") 

For  our  purposes,  yield  strength  can  be  considered  to  be  the 
stress  at  which  permanent  deformation  begins  to  occur,  and  tensil  strength 
the  stress  at  which  total  failure  occurs. 

To  ensure  a  margin  of  safety,  safety-factor  multiples  (i.e., 
design  factors)  are  customarily  applied  to  Instrument  and  machine  designs. 
That  Is,  the  design  Is  such  that  the  maximum  stress  in  the  device  will  be 
a  fraction  1/N  of  the  yield  strength  or  tensil  strength  of  the  construction 
material,  where  N  Is  the  safety  factor.  According  to  Falres^1^,  the 
following  are  traditional  values  of  N  for  steel  and  ductile  metals: 


11.  V.  M.  Fa ires,  Design  of  Machine  Elements  (MacMillan,  1965).  Table  1.1 


TENSI  L  STRESS  (psi  ) 


TENSIL  STRENGTH 


STRAIN  (  PERCENT) 


Figure  33.  Typical  stress-strain  curve 


Based  on  Based  on  • 

Yield  Strength  Tensil  Strength 

3-4  1.5-2 

6  3 

8  4 

10-15  5-7 

The  last  row,  for  shock,  would  be  the  one  applicable  here.  Thus,  the 
blades  should  be  constructed  of  a  metal  that  has  a  tensil  strength  5-7 
times  the  stress  values  computed  by  eqs.  (41)  or  (42),  or  shown  in  Figure 
32. 

Table  2  gives  yield  strength,  tensil  strength  and  some  other 
data  for  aluminum,  titanium,  several  strong  alloys  of  aluminum  and  titanium, 
and  for  a  strong  stainless  steel.  Consider  aluminum  with  a  tensil  strength 
of  13,000  psi.  This  should  be  5-7  times  the  maximum  stress  in  the  blades, 
and  from  Figure  32,  we  find  that  blades  of  unreasonable  thickness  would  be 
required.  Aluminum  alloys  are  considerably  stronger,  but  still  Inadequate. 
Titanium  alloys  and  stainless  steel,  with  tensil  strengths  on  the  order 
of  150,000  psi,  would  require  blade  thickness  of  about  3  mm,  including  a 
conventional  safety  margin,  to  withstand  collision  with  the  larger  drops. 


Kind  of  Load _ 

Dead  load 

Repeated,  unidirectional, 
gradual,  mild  shock 

Repeated,  reversed, 
gradual,  mild  shock 

Shock 


Sumnary  and  Conclusions 

Maximum  stresses,  caused  by  impact  with  large  water  drops,  in  the 
tall  portion  of  the  HOI  blade  design  shown  In  Figure  11  are  plotted  in 
Figure  32.  Maximum  stress  varies  with  the  square  of  the  ratio  of  water 
drop  radius  to  blade  thickness. 


Aluminum  and  aluminum  alloys  are  found  to  lack  sufficient  tensil 
strength  to  be  used  for  the  blades  (Table  2).  Titanium  alloys  and  strong 
stainless  steel,  which  have  tensil  strengths  on  the  order  of  150,000  psi, 
are  adequate.  Including  a  safety  margin  factor  of  about  5,  a  blade  thick¬ 
ness  of  about  3  mm  should  be  adequate. 
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APPENDIX  A 


GLOSSARY  OF  SYMBOLS 

water  drop  radius 

speed  of  sound  in  water 

defined  by  eq.  (29) 

particle  impaction  efficiency 

dynamic  pressure  force  on  a  torque  blade 

blade  height 

hydrometeor  spectrometer 
torque  blade  length 
water  drop  impact  pressure 
dynamic  pressure 

radial  coordinate  in  an  axially  symmetric 
system  with  the  symmetry  axis  in  the  x 
coordinate  direction 

inner  radius  of  the  annular  collection  surface 

outer  radius  of  the  annular  collection  surface 

radius  of  a  simple  disc  particle  collector, 
or  see  eq.  (28) 

defined  by  eq.  (4).  or  strain 

defined  by  eqs.  (12)  and  (27) 

torque 

nondimensional  torque  (eq.  (13)) 
flow  vector 


y(a,£l,S) 


maximum  impact  speed  of  a  water  drop  on  the  side 
of  a  torque  blade 

component  of  flow  normal  to  the  blade  surface 

radial  component  of  flow  velocity 

x  component  of  flow  velocity 

free  stream  air  speed 

cloud  water  content  meter 

streamwise  coordinate,  positive  in  the  down¬ 
stream  direction  with  origin  at  the  collection 
surface 

angle  between  the  torque  blade  and  the  radial 
direction 

angle  between  the  torque  blade  and  radial  axis 
which  defines  the  point  along  the  blade  at 
which  the  dynamic  pressure  is  zero 

angle  of  the  torque  blade  with  the  radial  axis 
at  r  *  r2 

torque  blade  thickness  (see  Fig.  31) 
air  density 
water  density 
defined  by  eq.  (16) 

critical  contact  angle  of  water  drop  during 
impact  on  a  rigid  surface  (see  Fig.  28) 

angular  rotation  speed  of  the  annulus 

equilibrium  angular  rotation  speed  of  the  annulus 

nondimensional  angular  rotation  speed  (eq.  (11)) 

equilibrium  nondimensional  angular  rotation  speed 


APPENDIX  B 


NEWTON-RAPHSON  CALCULATION  OF  EQUILIBRIUM  ROTATION  SPEED 

We  express  rotation  rate  in  terms  of  the  dimensionless  rotation 
speed  n  (eq.  (11)),  and  we  take  both  h  and  vf  to  be  constant.  Our  objec¬ 
tive  is  to  find  si  such  that  T  (eq.  (18))  is  zero.  We  begin  with  an 
initial  estimate  of  si,  si,,  and  after  i  -  1  Newton-Raphson  iterations  we 

X.  u  I 

calculate  the  i  +  ltn  estimate  of  a  from  the  equation 

ni+l  =  fii  '  ^./(dT/dsi) .  (1A) 

Iteration  continues  until  T  differs  from  zero  by  no  more  than  an  acceptable 
tolerance. 

Given  Q..,  we  can  calculate  from  eq.  (18),  and  it  remains 
to  derive  an  expression  for  dT /dsi.  This  can  be  shown  to  be 


_J _ 

sin2e 


dT 

dsi 


-2ctn<x  (  1 
o  \ 


ft  sine  etna  csca 
0  0 


2  f  4  4  4 

+H&  sin  e  2ctn  aQ  -  ctn  b  -  ctn  (<j> 


-sing 


2ctnao  _  ct.nB  _  ctnU  ±  b)  + 

sina  ’  sins  "  sin(ij>  +  b) 

0 


fcn 


tan  (a  /2) 

_ _ 0  _ 

tan  ( 6/2)  tan  (^-^—) 


where 


2 

da  tana  sin  a 

_ 0  _  _ 0 _ 0 _ 

dsi  3  2  /  ■  ~7 

n  sin  b\1  +  (2/sisinB) 


(3A) 
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